Gain-tunable optomechanical cooling in a laser cavity 
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We study the optical cooling of the resonator mirror in a cavity-optomechanical system that 
contains an optical gain medium. We find that the optical damping rate is vanishingly small for an 
incoherently pumped laser above threshold. In the presence of an additional external coherent drive 
however, the optical damping rate can be enhanced substantially with respect to that of a passive 
cavity. We show that the strength of the incoherent pump provides the means to tune the optical 
damping rate and the steady state phonon number. The system is found to undergo a transition 
from the weak optomechanical coupling regime to the strong optomechanical coupling regime as the 
strength of the incoherent pump is varied. 
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Cavity optomechanics is a rapidly growing field of re- 
search studying the couphng of mechanical degrees of 
freedom to modes of optical cavities [U [2]. One fasci- 
nating consequence of this interaction is the possibiHty 
to laser-cool the mechanical motion (e.g. [3HS])- Theory 
predicted ground state cooling to be possible [ZHSj, and 
recent experiments [101 E] reached down to less than 
one phonon. Most of the earlier work in this area has 
focused on optomechanics in a passive cavity setup with 
mechanical elements, with relatively little attention on 
non-standard setups including optical non-linearities |12j . 



In this Letter we address the following question: can 
a laser cool its own mirrors? More specifically, we con- 
sider a cavity-optomechanical system with an incoher- 
ently pumped optical gain medium, modeled as an en- 
semble of identical two-level systems that couples to the 
cavity photon field. We find that the answer is surpris- 
ingly not affirmative. While a laser can cool a second 
(passive) cavity, as studied extensively in earlier works, 
it is very inefficient in cooling its own mirrors. This is at- 
tributed to the saturating non-linearity that gives rise to 
a vanishing restoring force on cavity field fluctuations. In 
the presence of an additional coherent "seeding" signal 
however, we find that the damping rate can be signifi- 
cantly enhanced with respect to that of a passive cav- 
ity. We attribute this effect to a reduced effective cavity 
decay rate that can be tuned by the strength of the in- 
coherent pump. We subsequently discuss the minimum 
attainable phonon number and show that it can be low- 
ered from that of a passive cavity if the thermal contribu- 
tion is non- negligible. In the end we discuss the phonon 
spectrum and show a transition from the weak optome- 
chanical coupling regime to the strong optomechanical 
coupling regime as the incoherent pump strength is var- 
ied. 



We first consider the case without the seeding signal. 



The Hamiltonian can be written as 

H = ujr{l - ^)ia'a - n) + ujmb^b +'^y^ 

L 2 ^ — ' ■' 

j 

+ 3 + "^Z") + ^bath + ^pump, (1) 
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in which we have set fi = 1. a) ,a and , b are the creation 
and annihilation operators of the cavity photon field of 
frequency ujr and the mechanical eigenmode of frequency 
ujm , respectively. L is the length of the optical cavity, n 
is the average cavity photon number in the steady state, 
and X = xzPF(b-\- b^) in which the zero point fluctuation 
of the mechanical oscillator xzpF — (2tow„i)^2 . The first 
two terms in ([T]) describe the corresponding passive op- 
tomechanical system ^ ; the next two describe the active 
medium of transition frequency Ua and its coupling to the 
cavity field with strength g. We have employed the ro- 
tating wave approximation, and the subscript j indicates 
each two-level system in the ensemble. Below we use 
the collective operators P = , D = ct| whose 

expectation values give the macroscopic polarization and 
inversion of the medium, -ffpump describes the incoherent 
pump and the decay of the 2-level system it introduces, 
and TJbath represents all the other decay processes. 

A first insight into the modification of opto-mechanical 
interaction due to the atomic medium can be obtained by 
considering the optical damping rate Topt from Fermi's 
golden rule in the rate equation approach 7 : 

ropt('^in) = ^^[Snni^m) — Snn[—^m)\- (2) 

Here G = w^^zpf/^ is the optomechanical coupling 
strength, and Snniy) is the autocorrelation function of 
the cavity photon number, defined as the Fourier trans- 
form (F(w) = F{t)e''^^dt) of 

Snnit) = {a\t)a{t)a\Q)am - {a\t)a{t))\ (3) 

ropt(<^m) defined above applies to both the Hamiltonian 
(n]) and a passive system in which no gain medium is 
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present. The coupling between the cavity field and the 
gain medium brings the system into a self-organized op- 
tical oscillation in the presence of the incoherent pump, 
centered at the laser frequency lol. Consequently, the re- 
laxation dynamics of 5'„„(i) is modified with respect to 
that of a passive system and may display multi-mode or 
chaotic dynamics in general. 

Assuming the decay of the polarization (7j^) is much 
faster than that of inversion (7|| ) as in a Class A or B 
laser, we first adiabatically eliminate P in the rotating 
frame of the laser frequency lol (see Appendix|A]). By lin- 
earizing the operators about their classical steady-state 
values, i.e. a(t) = {a + M(t))e-*"^*, b{t) = D + SD{t), 
we obtain the following equation of motion 



Sa{t) = 



7j_ - iALa 



aSD{t)+Ta{t), 



(4) 



above the laser threshold and in the absence of the op- 
tomechanical coupling. Here |ap = n, A^a = — is 
the detuning from the atomic transition frequency, and 
J-a is the effective fluctuation force on Sa (see Appendix 
[A| ). The incoherent pump strength is quantified by the 
steady-state value of the inversion D in the absence of 
the cavity field, denoted by Dq. Above the laser thresh- 
old, i.e. Do > Dth = 2k/W, gain saturation clamps 
D to Dth P^HI5] . Here k is the cavity decay rate (of 
a(t)) and W = 2g'^"fj_/{^j_ + A'j^^) is the stimulated emis- 
sion rate. The laser frequency is determined by the line- 
pulling formula Air + '^Aia = 0, setting a relationship 
between the atomic detuning Ai^a and cavity detuning 

Alt = UJl — {jJr- 

A crucial feature of Eq. Q is the vanishing of the 
complex restoring force on dd{t), i.e. the term propor- 
tional to 5a{t) on the right hand side. This can be 
attributed to the saturating non-linearity of the laser. 
The cavity decay as well as the detunings precisely can- 
cel by virtue of the saturating gain and the line-pulling 
formula. The photon autocorrelation function is conve- 
niently calculated by expressing Snn(t) — {6n{t)6n(0)) in 
terms of the photon number fluctuation operator 6n(t) ~ 
a6a^ (t) +a*6a{t) and solving the coupled Langevin equa- 
tions for the pair {Sn,SD) (see Appendix LaI): 



+ 7^(1+0 



2- WNg n. 



(5) 



Here Ng is the total number of gain atoms, ^ = Dq/D ■ 
1 = n/usat is the dimensionless saturation factor, risat - 
7||/2W is the saturation photon number, and uj± : 



;Tl| 



(1 + ±^(1 + 0^ 



8C 



7|| 



are the complex re- 



laxation frequencies of the laser. Eq. ^ shows that 
the optical damping rate Fopt vanishes because 5'„„(w) 
is symmetric with respect to w = 0. This result is valid 
for ^ < 1, 7|| /7_L ^ 1 (see Appendix [A|) and can be traced 
back to the vanishing of the complex restoring force on 
Sait). 



Within the rate equation approach, the steady-state 
phonon number can be written as fim = + 
G'^ Snn[~u]ra) m for Fopt — (scc Appendix p3|) , where 
F,„ is the mechanical damping rate and nj-^ is the phonon 
number in thermal equilibrium with the mechanical bath 
alone. The result shows that the mechanical oscilla- 
tor acts as a spectrometer for the intracavity radiation 
pressure noise generated by the active medium. It can 
be shown from Eq. ([5| that 5'„„(— w™) is positive and 
Uyn > n^i the optomechanical coupling increases the 
effective temperature of the mechanical motion, due to 
the additional noise introduced by the photonic and the 
atomic medium. The same conclusion can be drawn from 
the more rigorous approach of integrating the phonon 
spectrum S'bt6(w) = (S^(t)fe(0))e'"* (see Appendix 

0- . . . 

The situation is vastly different when a seeding signal 

of frequency is fed into the cavity coherently, in ad- 
dition to the incoherent pump on the atomic medium. 
The seeding is described by a time-dependent Hamilto- 



-\- h.c), where fig denotes the 



nian Hg 

strength of the seeding signal. Note that wz, of the exter- 
nal laser drive can be chosen differently from the above- 
threshold laser frequency determined by the line-pulling 
formula. The system Hamiltonian in the rotating frame 
of frequency ujl becomes 

H = -[ALr + G{P + b)]{a^a - n) + Lj,J^b + y D 
+ ^g{a''aj +a+a) + Qs{a'^ +a). (6) 



The average cavity photon number n is now fig 
dependent as well: 



9'D 



^^A, 



'«)(7-L - iALa) 



iA 



Lr 



(7) 



where the steady-state inversion D = Do/(l -|- ^) as be- 
fore. Eq. Q can display a bistable behavior [H], but we 
will focus here to the resolved side band limit (k <C | Airl) 
with a small C 1; where a single steady-state solution 
exists. 

The photon field fluctuation now follows 



Sa{t) — (iA^r — K)6a(t) 



g^aSD{t) 



1± 



iA 



La 



+ Mt), (8) 



where the effective detuning A^r and cavity decay k are 
defined as 



ALr = A 



Lr 



WD 

2^ 



Aia, k 



WD 



(9) 



Note that D has the same sign as Dq, and for Do > 
this leads to an effective cavity decay that is reduced 
from its intrinsic value. This tunability of the effective 
complex frequency of photon fluctuations stems from the 
fact that D is no more clamped at Dth, as in the case 
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above the laser threshold in the absence of seeding. In- 
stead, it's now a function of the incoherent pump Dq and 
the coherent drive amplitude jri^j. Eq. ^ leads to the 
following form of Snn (see Appendix [A|): 

(10) 



where we have defined A(a;) 
e)2]/|e(w)|2 and Q{uj) = j^[{-iuj + k)^ 
J^^{^+0]-'2nW^D[ALaALr-l±ifi-^Lo)] 



The previously 

studied cases can be recovered in various limits of this 
expression. For the case without seeding k = A^r = 0, 
Q{ui) becomes —iuj^^(uj — a;+)(w — aj_) and we recover 
S'„„(w) in Isj) . T he passive case is recovered by setting 
f; = in Eq. ( |10[ ), reproducing the well-known expression 
5„„(w) = 2Kkj\{ijj + AlvY + K^] [m. 

Snn{^) given by (10) is asymmetric around w = and 
leads to a nonvanishing Fopt . For a moderate cavity pho- 
ton number, Qiuj) is dominated by the first term and 

A(w) W A(-W) « + [UJ- ALrf][k'^ + {UJ + KLr)\ 

indicating that the relaxation frequencies of the seeded 
laser are roughly iA^r — As a consequence, Snn{^) 
peaks near iA^r with a width given approximately by 
k (see Fig. [ijb)). These peaks are higher than that in 
the passive case due to the atomic noise, resulting in the 
optical damping rate 



opt 



-8G-^ nujAi^rfi 



[k2 + Ai,)2][K2 + + Ai,)2] ' 



(11) 



In the resolved side-band limit (of the active case, k <C 
|Aij.|), the maximal optical damping is therefore 



To 



pt, max 



/ k 



-\ -1 



(12) 



This results in a gain-enhancement of the damping rate 
by about n/k, when compared to the corresponding pas- 



sive system with the same mean number of cavity pho- 
tons (see Fig. [T]) Jjj. The maximal damping rate (12) 
requires a negative detuning A^^ = ^'^in if the effective 
cavity decay k is positive, similar to a passive system. 
If the incoherent pump is strong enough that k becomes 
negative, a positive detuning A^^ = is required in- 
stead, otherwise Fopt becomes negative, resulting in op- 
tical heating. We note that the mean number of cavity 
photons mostly is set by the strength of the coherent 
drive in the resolved sideband limit (see Eq. ([7|), as in 
a passive system. Yet the damping rate can be greatly 
enhanced due to a strongly modified dynamics of the pho- 
ton number fluctuations. 

We note that the enhancement of Fopt does not nec- 
essarily lead to a reduced temperature of the effec- 
tive optical bath, which can be represented by n^* — 
{Snn{u!m) / Snn{-i^7n) - 1)"^ (sce Appendix |b]) . The to- 
tal phonon number is given by h„i — (Foptf^m' + 
rm"-m)/(ropt + r„), and n^'(< ^m) is the mini- 
mum phonon number, attained in the limit Fopt 2> 
Tm, rm^m/'^m'*- l^is limiting regime the addition of 
gain leads to a higher effective temperature of the me- 
chanical motion, since n$^* is increased due to the en- 
hancement of both Snn{^v 
noise: 



,), Snn{-^m) by the atomic 



= opt 



WNg + 2{k - k) 
4k 



(13) 



However in the regime F^ ^ Fopt ^ Tm^-m/^'m' lhat is 
common in passive systems, either due to a small n^*, 
a relatively high thermal temperature, or a small G, the 
thermal contribution Fmn^/Fopt in cannot be ne- 
glected in general. In this parameter regime we find that 
the cooling of the mechanical motion can be more ef- 
fective with the addition of the gain, due to the stronger 
suppression of the thermal contribution as a consequence 
of an enhanced Fopt- 

One such example is given in Fig. |2j in which fi^* is 
fixed and the weight of the optical and thermal contribu- 
tions in fim is adjusted by considering different values of 
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FIG. 1. (Color online) Optical damping rate (a) and pho- 
ton autocorrelation function (b) in a seeded laser (solid line) 
and the corresponding passive system (dashed line). Param- 
eters in this typical example are: A_ta = — Al,. = 1 GHz, 
7x = 10 GHz, 7|| = K = 100 MHz, g = 1 MHz, Do = O.SiVg = 
1.2Dth, and n = 10^ The resulting ^ = 0.396, k/k = 7.12, 
and Fopt, max is 6.94 times its value in the corresponding pas- 
sive system. 



FIG. 2. (Color online) Steady-state phonon number versus 
the optomechanical coupling in a seeded laser (solid line) and 
the corresponding passive system (dashed line). The dotted 
and dash-dotted lines indicate the contribution from the opti- 
cal bath and the mechanical bath in the seeded laser, respec- 
tively. In the passive system n™ comes almost all from the 
latter, nj^ = 10^, LOm = —^Lr, Tm/K = 2 X 10~*, and the 
other parameters are the same as in Fig. [l] 
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FIG. 3. (Color online) (a) Steady-state phonon number ver- 
sus incoherent pump strength in a seeded laser (solid line). 
The dotted and dash-dotted lines indicate the contribution 
from the optical bath and the mechanical bath, respectively. 



The circle and squares show the approximation ( 14 1, ( 15 1 and 



the result from integrating the phonon spectrum, respectively. 
The dashed line indicates the phonon number in the corre- 



sponding passive system. G/k, = 3 x 10 



= 10'', and the 



other parameters are the same as in Fig. |2] k becomes zero 
at Do/Dth = (1 -I- ^) — 1.396. Inset: Optical damping rate 
in the seeded laser [18]. (b) Minimized phonon number at 
different mechanical frequencies via the optimization of Do- 
Air = —idrn IS adjusted accordingly. Dots show the fitting 



using Eq. ( 15 1 multiplied by [1 + 6(k/cj„ 



G. For G/k < 3 x 10"'* the seeded laser studied in Fig. [l] 
has a lower phonon number, which has a minimum close 
to fi°P* ~ 9. We note that the upper limit of this range 
scales linearly with {n^/n)^ . 

For a fixed G and n, the incoherent pump strength Dq 
provides a means to minimize n„i- If the optical or ther- 
mal contribution in rim is dominant for all possible values 
of Dq, the optimization is straightforward and achieved 
with the lowest nj^' or the largest enhancement of Fopt , 
respectively. When this is not the case, the optimiza- 
tion is achieved when nj^' roughly equals the thermal 
contribution, i.e. 



G^n{WNg + 2k) 

mii'm 



(14) 



at which the phonon number is approximately their geo- 
metric mean 



[WNg + 2K)Tmnl 
2G^n 



(15) 



We have taken a;„j = —Aj^^ and neglected the weak w„ 
dependence in (14) and (15), which is only second order 
in k/ujm as illustrated in Fig. |3i 



This implies that the 
optimization of Dq is achieved almost simultaneously for 
all the mechanical frequencies in the resolved sideband 
limit. 

rim calculated by the rate equation agrees well with 
results obtained by integrating the phonon spectrum, i. e. 
rim — 5^ J duj {b\u})b{—u!)) , as shown in Fig. pla). We 
note that aG ~ 0.0 Iw^ in this example, and it is more 
than one order of magnitude smaller than kVm for a small 




0.4 0.8 1.2 
Incoherent pump Dq/Dj^, 

FIG. 4. (Color online) Phonon spectrum logj^Q(S(,t6) for the 
seeded laser in Fig. [3[a), showing optomechanical normal 
mode splitting as a function of the incoherent pump strength. 
The system is in the weak optomechanical coupling regime in 
the absence of the gain (i.e. Do = 0). 



Dq. This places the system in the weak optomechanical 
coupling regime, and in the phonon spectrum there is 
only one peak near uj — —ujm (see Fig. |4]) . The reduction 
of fim as Dq increases from results in the lowering and 
a slight broadening of this peak [12] . As Dq is further 
increased, k ^ and the system undergoes a transition 
to the strong optomechanical coupling regime, since kF,„ 
becomes comparable and even smaller than dG, which is 
kept fixed in this process. Consequently, the peak at a; = 
—ujm in the phonon spectrum splits into two, resulting in 
hybrid resonances separated by roughly 2dG [7: ; their 
enhancement with Dq reflects the increase of fim in this 
regime. 

In summary, we have discussed the optical cooling of a 
cavity mirror in the presence of an incoherently pumped 
atomic medium. In the absence of an external laser drive, 
the coherent radiation in the cavity generated by the in- 
coherent pump above threshold results in a vanishingly 
small optical damping rate. This always results in the 
heating of the cavity mechanical oscillator due to the ad- 
ditional photon noise generated by the atomic medium. 
In the presence of an additional external coherent seed- 
ing, we show that the cooling rate and the minimum 
attainable phonon number can be tuned through the 
strength of the incoherent pump, which can be enhanced 
and lowered respectively from their passive cooling val- 
ues. In addition, we showed that the strong optome- 
chanical coupling regime can be reached by increasing 
the incoherent pump strength, even at a fixed number of 
cavity photons. 

We acknowledge helpful discussions with Girish Agar- 
wal. L. G. and H. E. T. acknowledge NSF for Grant 
MIRTHE EEC-0540832. F. M. acknowledges DARPA 
ORCHID, an ERG Starting Grant, DFG Emmy-Noether, 
and the ITN cQOM. 



Appendix A: Equations of motion 

The optical field in the cavity couples to the mechan- 
ical motion of its mirror and the polarization of the op- 
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tically active medium with strengths G = ^^^cjr and 
g = — /i^ojr/eoV^, respectively: 



a = {—iujr — K)a — igP + iGib^ + b)a + J^^ 



(Al) 



K is the cavity decay rate, /i is the dipole strength, cq 
is the permittivity of the vacuum and V is the effective 
volume. J^'s here and below are the fluctuating Langevin 
forces associated with the corresponding decay channels. 
In the rate equation approach the mechanical backaction 
on the dynamics of the rest of the system is neglected, 
which we will restore in the discussion of the phonon 
spectrum. 

The polarization P and inversion D follow 

P = {-iuja --/A_)P + igDa + P^^ , (A2) 

D = 7,1 {Do -D) + 2ig{a^P - P"^ a) + P^^^ . (A3) 

The incoherent pump strength is characterized by the 
parameter Dq which is the steady-state value of the in- 
version created in the absence of a cavity field, and 7|| 
are the decay rate of the polarization and inversion, re- 
spectively. Assuming 7|[ ^ as in a Class A or B laser, 

we adiabatically eliminate P in the rotating frame of the 
laser frequency lul [T^ : 



igaD + T- 



7-L 



7-L 



(A4) 



The linearization of the operators a, a\ D about their 
steady-state value a, a*, D in the same rotating frame 
leads to 



Sa = (iA^r — K)Sa + 



g^[DSa + a6D] 



Pa, (A5) 



7_L - i^La 

where the effective fluctuation force on da is given by 



T = T — 



7-L 



*A,,„-^-^-- 



(A6) 



Pa satisfies {Pa{t)Pi{t')) = 2?,,t<5(t-t'), (-^1 W-^a(i')> = 
^ata^(^ ~ ^'); whcrc the diffusion coefficients P^at = 
^Lt t and = -Dlfa+T^lfa contain a blackbody 
contribution + — 2K(n^ -1-1), = 2Kri^ and a con- 
tribution due to collective atomic dissipation processes 



SE 


w 


a'^ a 


~ Y 


SE 


w 




~ T 



D) 



27- 



(Do-D) 



(A7) 
(A8) 



is the number of blackbody photons inside the cavity, 
vanishingly small at optical frequencies studied here and 
Ng is the total number of gain atoms. 

Eq. ( |A5[ ) is further simplified to (4) in the main text 
by using the values of D and a. Similarly, we find 



where n is the average cavity photon number and W = 
25^7.1/(71 + A|^) is the stimulated emission rate. The 
fluctuation force P~^^^ satisfies {IF-y^^F^^^) = 'D^^^~f^^5{t — t'), 

(PaP^,) ^Va^^^sit^f), (j-t^^,,) = 2?„t^„5(<-i') 

where 



V. 



7||7|| 



2711 N, 



V. 



«7|| 



N„ 



l\\g^Da 
(7± - *Ai„.)2 
7||, g^ZJa 
(7_L + iALa)'^ 



Do 
N, 



9 

Do 



(AlO) 
(All) 
(A12) 



The equations of motion for Sn{t) = a6a^{t) + a*6a{t) 
and SD(t) in the frequency domain can then be written 



M- 



dPl{uj)+a*PaiL0) 



7|| 



where 



M2 



( -e7||/2 
\2WD -iw + 7||(l + e) 



(A13) 



(AM) 



Here ^ = Do/D — 1 = n/usat is the dimcnsionless sat- 
uration factor, which is a convenient parameter used in 
the main text to measure the number of photons gener- 
ated, and Us at = ■^ \\/2W is the saturation photon num- 
ber. Solving (A13) we find 



6n{u!) 



zo^- 711 (1 + 
(cj — — a;_) 

2{u! — u;+)(w — UJ-) 



7|| ' 



(A15) 



which lead to 



5'riTi('^) 



0.2+72(1+02 

2_,„2U,„2_,.,2N " (^-t + 2?ata) 



(a;2 -t^2)(^2 _^2) 

lie 



4(a;2 -a;2)(cj2 _^2) 



2 ^ -^711 7|| 



(1+0^7|5 

(up- — ijj\)(yp- — Lo 

w^7||a 
(cj2 _a;2)(w2 _^2) 



2 \ -R'*^[^a7|| 



Im[I?a 



(A16) 



using the Wiener-Khinchin theorem Snn{^) ~ 
{Sn{uj)Sn(—Lu)) . We have set a = a* by choosing a 
proper t = 0. uj± are the complex relaxation frequencies 
of the laser as defined in the main text. 



We first note that the first terms in Eq. (AI6) is much 
larger than the rest for ^ < 1, 7||/7j_ ^ 1, with which 
Eq. (5) in the main text is derived. Furthermore, we 
have neglected the blackbody contribution to the noise, 
SD = -(7, +2Wn)5D~2WD[dSa'' +a*Sa]+P^ (A9) i.e. taking V^ia + 'Daaf = WNg + 2K{2n^ + 1) « WNg 
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in the first term. It is symmetric with respect to w = 0, 
and so are the second and third terms in Eq. (A16). The 
asymmetry of iS'„„(aj) only arises from the last term and 
is vanishingly small. This derivation shows that only 
leads to a small effect, and we have thus neglected it in 
the rest of the discussions in the main text. 

The main approximation we have employed in deriv- 
ing Eq. (5) in the mai n te xt is the adiabatic elimina- 
tion of P{t) (see Eq. (A4|). This is typically a good 



approximation for a Class A or B laser. To demonstrate 
its validity in the calculation of Fopt, below we compare 
our results to the solution of the full linearized equa- 
tions of motion for the photonic and atomic degrees of 
freedom, i.e. keeping the fluctuation SP{t) defined by 
p{t) = (P + (5P(t))e~*"^*. SP{t) evolves according to 

6P = {iALa - 1±)SP + ig{D6a + a5D) + P^^ . (A17) 

By solving it together with its hermitian conjugate and 
Sa, Sa\5D, we again derive an expression for S'„„(a;). 
Due to its complexity, here we only compare it with 
Eq. (5) numerically. As Fig. [s] shows, the asymmetry 
of Snn{ijj) about u — Q arising from the fast dynamics of 
5P{t) is extremely small, so is the optical damping rate 
Fopt- 



10^ 



CO 

10-^ 



-1 -0.5 0.5 1 
a)(GHz) 

FIG. 5. Photon autocorrelation function calculated with 
the polarization P{t) adiabatically eliminated (solid) and lin- 
earized (dashed). They agree well, with the asymmetry of the 
latter (i.e. |5'„„(a;) — S'„„(— aj)|) shown by the dash-dotted 
line. Parameters used are: 7x = 10 GHz, 7|| = k = 100 MHz, 




Ai,„ = 10 GHz, g = l MHz, n = 10^ and Dt 



: 0.17V,, 



The vanishing of the complex frequency of Sa^{t) 
Eq. (4) in the main text decouples the dynamics of 6a 
and Sa^ for the unseeded case. The situation is different 
with a non-vanishing seeding signal, where the phase is 
enslaved by the phase of the coherent driving field. We 
therefore introduce du(t) = a*da{t) — ada^{t) along with 



6fi{t) and 5D{t). Their equations of motion with the 
Hamiltonian (6) in the main text can be written as 



6u{uj) 

where 

~iuj + k 

M3 = I -iALr 

2WD 



-iuj + k 




-Wn 



(A18) 



(A19) 



Their solution leads to the autocorrelation function (10) 
in the main text. To recover Eq. (A13) in the unseeded 



case, we note that Sii is decoupled from Sri and SD since 
k — A^r = 0. 

Next we discuss the calculation of the phonon spec- 
trum 5*616 = (6^(w)&(— oj)) with the linearized equations 
of motion. The dynamical equation for the mechanical 
motion is 



{-iUJn 



)b + iG{a'a~n)+J^r„ 



(A20) 



where ujm and F„i are the frequency and the intrinsic 
damping rate of the mechanical oscillator. Solving it 
together with 6a{t), 6a^t) (or Sh{t)) and SD{t) in the 
frequency domain, we find in general 

(A21) 

For an unseeded laser, we find that Xfet = 0, Xm(w) = 
[i{ujm Fm/2]^^ are exactly the same as in the case 

of a standalone mechanical oscillator, as a consequence 
of the vanishing of the complex restoring force on 6a(t). 
The latter also leads to Xa(w) — Xa^i^), given by 



aG[^ + »7||(l + /3)] , , 

Xa W) = rXm W . 



(A22) 



The resulting phonon number = ^ /T^^ Si,-tf,{uj)du} is 
then the sum of the unaltered thermal contribution from 
the mechanical bath. 



-T _ 



1 

2^ 



|Xm(-w)pX'bt6dw = 



(A23) 



and the extra contribution form the optical noise 
^IZ. du;\xai-u;mWNg + 2^i2n^ + l)]>0. 

In a seeded laser the response functions are more com- 
plicated, and here we only give Xa(w) as an example 



Xa(w) 



iaG'x;(-cj)xrn(t^) 



1 -H4nG2cJ™ALrXr(-w)Xr(w)Xm(--^)Xm(t^) + Zjfk 



2inWDg'^ 



1+«)7|| 



X;(-") _|_ Xr-i'-') 



(A24) 



where Xr{^) = [iiA^r ^ oj) + k] ^ . The transition of the system from the weak optomechanical coupling regime 



7 



to the strong optomechanical couphng regime can be 
analyzed from the denominator of Eq. (A24). We first 
neglect the last term, which is absent in the passive 
system. The hybrid resonances appear when the real 
part of the denominator becomes zero, which occurs at 
cj « —uJm ± aG and requires that aG is comparable or 
larger than kT„i . The last term and the minute difference 
between A^r and A^,. slightly change the separation of 
the hybrid resonances, as well as making them asymmet- 
ric about W — —LOra- 



Appendix B: Rate equation model for mechanical 
oscillation 



In the main text we have used a rate equation ap- 
proach, where the rates are calculated in the quantum 
noise picture as in Ref. [7]. This approach is a good ap- 
proximation if the relaxation of the effective optical bath 
(e.g. K in a seeded laser) is much faster than its energy 
exchange rate with the mechanical oscillator (i.e. Fopt), 
with which we can treat the former as a Markovian bath. 
This condition is satisfied for Dq/Dih ^ 1 in Fig. 3(a), 
and its prediction is qualitatively correct even beyond 
this range as we find. 

Denoting p„ „ as the occupation probability in the 
state with n phonons and F„^„-|-i(a;m) as the transition 
rate from the state with n phonons to the state with n± 1 
phonons, the detailed balance in equilibrium requires 

(Wm) + P'n—l.n—l^'n—l.n 

M = 0, (Bl) 

i.e. the transition from the state with n phonons to the 
state with (rt — 1) phonons are balanced by the reversed 
transition. Using the system-bath theory or Fermi's 



golden rule, we find that r„.„_i and r„_i^„ are propor- 
tional to n, and we denote Fi o = F^ and Fo,i = F^. 

Assuming the equilibrium satisfies a thermal distribu- 
tion, one can easily show that the average phonon number 



Ft 



F 



(B2) 



If the transitions are only caused by the intrinsic me- 
chanical damping, the average phonon number is just 
the thermal phonon number n^. Similarly, we can de- 
fine a nj^* if the transitions are only caused by the optical 
damping, and its effective temperature is give by 



2opt — 



hun 



[ln«r + 1) - ln«r)] , (B3) 



where fc^ is the Boltzmann constant. With both damp- 
ing processes present, fim can be rewritten as a weighted 
average of these two phonon numbers 



i opt I I'm 



opt 



(B4) 



n^* becomes ill defined if Fopt vanishes, as in an un- 
seeded laser. In this case n,„ can be expressed as 



popt 



G'^ Snni^^m) 
f 

^ m. 



(B5) 



It has the same structure as the result from integrat- 
ing the phonon spectrum S^f),, i.e. the unaltered ther- 
mal photon number plus a contribution from the photon 
noises. We note that here the relaxation rate of the opti- 
cal bath is not k(« 0) but rather — Im[c<;±](7^ 0), and the 
rate equation still applies. 
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